General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 






f-; Mu ¥f 3/?a7* 

(U<k. : — 


March 22, 1983 


SRI International 
Final Report 


THEORETICAL STUDIES OF PHOTODISSOCIATION OF SMALL 
MOLECULES OF ASTROPHYSICAL IMPORTANCE 


By: Roberta P. Saxon 


Contract NASW-3373 
SRI Project No. 1430 
MP 83-063 


Prepared for: 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Earth and Planetary Exploration Division 
Washington, DC 20546 


Attn: Dr. Henry Brinton 




Approved: 



Jonald C. LorentTs/ Director 
Molecular Physics laboratory 


G. R. Abrahamson 

Vice President ‘ 

Physical Sciences Division 

(NASA—CB— 1702 15) THEORETICAL STUDIES OF N83-23262 

PHOIODISSOCIATI ON OF SMALL MOLECULES OF 
ASTROPHYSICAL IMPORTANCE Final Report (SRI 
International Corp - f Menlo Park, Calif . \ Uaclas 

13 p HC A02/MF A0 1 jCSCL 03B G3/9 1 10054 

333 Ravenswood Ave. . Menlo Park, CA 94025 1 

{415)859-6200 . TWX; 910-373-2046 . Telex: 334 486 r 



CONTENTS 


*? 


BACKGROUND . 1 

RESEARCH ACCOMPLISHMENTS * 3 

Photodissociation of NH 2 ...... . . » • < ...3 

Calculations on C 2 H. . . ......6 

REFERENCES 


APPENDIX: 

Theoretical Study of NH 2 : Potential curves, transition moments and 

photodissociation cross sections 


i 


BACKGROUND 


It was suggested many years ago that the radicals and ions observed in 
comets result from photodissociation and photoionization of molecules. 
According to current models/ a comet is composed chiefly of a large, solid 
nucleus of frozen gases (parent molecules) such as 1^0, HCN, and N^. It is 
believed comets were formed at the same time and in the same region of space 
as the major planets and that their chemical composition may be the same as 
that of the early solar system. As the comet nears the sun, the surface heats 
up, liberating the frozen gases as well as dust particles. Solar radiation 
photodissociates the parent molecules into fragments that can be observed by 
resonance fluorescence. 

Both polyatomic molecules, known to be present in the interstellar 
medium, and cometary radicals have been observed. Using laboratory photo- 
dissociation data and computer models, astronomers are attempting to identify 
the parent molecules that account for all observed radicals and ions. In many 
cases, they have been very successful. The water molecule is the likely 
precursor for the H, OH, and O^D) species that have been observed. Simi- 
larly, evidence suggests that HCN and CH^CH are the precursors for the CN 
radical.^ 

For some fragments, however, the precursor molecule is not as clear. One 
of these is the C ^ radical, which is observed by the Swan band emission 

O "3 

d^ll^ •* a II,,. This band system is one of the strongest observed in comets. 

One likely candidate for the C£ parent species is the acetylene molecule, the 

photodissociation of which has been studied experimentally in the vacuum 
6 7 

ultraviolet. * Absorption of a single photon produces excited acetylene that 
may dissociate to the C 2 H radical and an H atom. For wavelengths less than 
1200 A, the quantum yield for the latter process becomes significant. No evi- 
dence for production of the C 2 radical by absorption of a single photon by 
C 2 H 2 has been found. Using the ArF laser at 1930 A, researchers have 
observed the C 2 (d ILg ) emission, which is the result of a three-photon process 
in C 2 H 2 that may proceed through a C 2 H intermediate. 
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It has been suggested^ that the C 2 Swan band emission observed in comets 
Is a result of the secondary photolysis of the C 2 li radical produced by photo- 
dissociation of acetylene* This hypothesis is supported by the observation of 

O 

the C 2 H radical in interstellar space. 

The parent molecule of the NH radical is also ambiguous. The most likely 
parent for NH is ammonia, which has been observed in interstellar space but 
not in comets. The photodissociation of NH3 has been studied in the labora- 

Q 

tory in the vacuum ultraviolet. The primary product of dissociation is known 
to be ground state + H. Fluorescence from excited state NH/^CA^A^) has 
been observed^ for incident wavelengths between ~ 1100 and 1300 A. Triplet 
emission of NH(A^II X^2) has been observed^ in comets, although only singlet 
NH is found as a direct result of NH3 photodissociation. Dissociation of 
ammonia to give triplet NH + ^ is a spin-forbidden process, and dissociation 
of NH3 to give triplet NH + H + H requires a great deal of energy. 

The suggestion that secondary photolysis of the N^ radical is related to 
the observation in comets of triplet NH originally motivated us to study the 
photodissociation of NI^. However, calculations on N^ alone will not 
unambiguously determine the parent molecule of the NH triplet state because of 
the following facts. The lower state of the observed NH triplet emission, the 
X E state, is the ground state of the system. Given a lack of competing 
processes, NH excited singlets will eventually radiate to the triplet ground 
state even though single-triplet transitions are dipole "forbidden." The 
observations do not indicate in which state the NH molecule was initially 
formed. 

We therefore broadened our goal to exploring NH*, photodissociation pro- 
cesses leading to not just triplet states, but also to singlet and triplet 
states of NH. Resulting photodissociation rates as a function of incident 
wavelength may then be incorporated into computer models of comet chemistry 
and physics. 

In this final report, we describe our calculations on NH2 and C 2 H 
performed under NASA contract NASW-3373. 
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RESEARCH ACCOMPLISHMENTS 


Calculations wore performed first for tho NH 2 system because it is 
computationally easier to treat than C 2 H since it contains two hydrogen atoms 
and only one first-row atom. The theoretical approach was carefully tested. 
Results were used in the calculations on C 2 H. 

Photodissociacion of NH 2 

Technical details of the calculations on NH2 have been published in the 
Journal of Chemical Physics; a copy of the paper 11 is attached as the 
Appendix. In this section we will briefly summarize the important steps in 
the calculation and the results of interest for cometary modeling. 

Spectroscopic information is available only for the ground X(^B^) and 

first excited A( 2 Aj) states of NH 2 ; therefore, the first Important task was to 

establish the correlation diagram for NH 2 dissociating to NH + H* Assuming 

the Franck-Condon principle, excited state energies ware determined near the 

equilibrium geometry of the ground state. If we consider only doublet states, 

which are dipole connected to the ground state, the first three states of NH 2 

are valence in character and the next two doublets are the lowest 3s Rydberg 

states in Bj and symmetry. Vte limited our consideration to these states, 

1 1 a* 3 

which include photodissociation channels to the a A, b £ , and A II states 
of NH. 

Potential energy curves for these five states as a function of one NH 
bond length at fixed remaining NH bond length and bond angle were performed 
using the density-matrix-directed second-order Multi-Configuration 
Self-Consistent Field (MCSCF) procedure of Lengsfield. A separate calculation 
was performed for each state. Particular care was taken to define a single 
calculation that describes equivalently pure valence, pure Rydberg, and mixed 
valence-Rydberg states. Extensive testing of this procedure was performed. 
Transition dipole moments were evaluated from the MCSCF wave functions. 

Using the calculated potential curves and transition moments, we calcu- 
lated photodissociation cross sections as a function of incident wavelength 
for transitions from the ground state, treating NH 2 as a pseudodiatora HN-H. 
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If we assume that each photon absorbed leads to dissociation, the cross 
section for photodissociation is given by 

o(v) - i F v |<x v |li|x k >| 2 ( 1 ) 

V 


where v is the frequency of the transition, X y is the bound state radial 
wavefunction, X k is the dissociating state radial wavefunction, p is the 
transition dipole moment, and F v is the fraction of molecules initially in 
vibrational state v. We do not consider rotation in this work. Although a 
Boltzmann average over initial vibrational levels is Indicated in Eq. (1), in 
this system at 300 K, only the ground vibrational level is significantly 
populated. 


The spectroscopic parameters of the pseudodiatomic ground state are 
reasonably similar to the parameters of a three-body system due to the small 
mass of the hydrogen atom. For example, the calculated w e value for the 
ground state pseudodiatom is 3507 cm which may be compared with the oalcu* 

O 

la ted symmetric stretch and asymmetric stretch of NI^CX ^B^), 3307 and 
3627 cm -1 , respectively. 12 Therefore, averaging over the bound state 
wavefunction in Eq. (L) is at least qualitatively reasonable. 


Photodissociation cross sections for transitions from the ground 

vibrational and rotational level of the 1 4 A" state to the 2 A’ and 3 A ? 

states are given in Table 1 for wavelengths in the center of the grid used by 

Huebner and Carpenter^ and are plotted in Figures 5 and 6 of the Appendix. 

Absorption to the 2 2 A' state results in a two-peaked cross section with peaks 

—21 2 

centered near 1710 and 1896 A with maximum vales of 2.7 x 10 cm and 
4.3 x i(T 21 cm 2 and halfwidths of 132 A and 68 A, respectively. The structure 
in this cross section is a consequence of the zero in the transition dipole 
moment function (due to symmetry) at the symmetric C 2 V geometry. The very 
small nmgnltude ot the cross section reflects the fact thar. the 2 *A’ state is 
bound with respect to dissociation in the center of the Franck-Condon 
region. Thus, the cross section arises from only the very edge of the 
Franck-Condon region. The cross section for transitions to the 3 A* state 
has its main peak at a value of 2.9 x 10“ 17 cm 2 with a halfwidth of 9 A at 
1622 A. The structure in this cross section is due to the transition moment 
crossing zero. 
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Table 1 


nh 2 PHOTODISSOCIATION CROSS SECTIONS* (cm 2 ) 


iV 

- 2 V 

>.d) 

a 

1526-. 5 

2.617F.-23 

1550.0 

1 .25ZE-22 

1574.5 

3.623E-22 

1600.0 

7 » 991E-22 

1626.0 

1 » 429E-21 

1653.0 

2.150E-21 

1631 .0 

2 ♦ 694E-21 

1709.5 

2 . 712E-21 

1731.5 

2 . 249E-21 

1746.5 

1 . 726E-21 

1762.0 

1.096E-21 

1773. 0 

4.S95E-22 

1794.0 

8 » 426E-23 

1310.0 

2 . 775E-23 

1326.5 

4.340E-22 

1343.5 

1 . 326E-21 

1860.5 

2.473E-21 

1373.0 

3.532E-21 

1396.0 

4 ’* 292E-21 

1914.0 

3.905E-21 

1932.5 

1 .S45E-21 

1951.5 

5.031E-23 

1 970 . 5 

2 . 753E-25 


1 2 4 

« „2 / 
. - 3 A 

\ ch 

a 

1262.0 

1. SUE-21 

1270.0 

2 . 365E-21 

1278.0 

2.973E-21 

1236.0 

3 ♦ 6I2E-21 

1294.5 

4.293E-21 

1303.0 

4 ♦ 891E-21 

1311.5 

5 . 330E-21 

1320.0 

5 ♦ 557E-21 

1328.5 

5.439E-21 

1337.5 

4 ♦ 921 E— 21 

1346.5 

4 . 004E-21 

1 355 ♦ 5 

2 . 731E-21 

1360.0 

2 . 04 IE— 21 

1 374 . 5 

2 . 073E-22 

1334.0 

1 ♦ 660E-22 

1403 . 5 

7. 536 E- 21 

1418.0 

2 » 451E-20 

1433.5 

7.234E-20 

1459.5 

1 .404E-19 

1481.0 

2 • 26 9 E- 19 

1503.5 

2.453E-19 

1526.5 

1 .361E-19 

1550.0 

1 .675E-23 

1574.5 

5 . 339E-1 9 

1600.0 

3 * 587E-1S 

1626.0 

1 ♦ 854E-17 

1653.0 

1 ♦ 269E-19 


Calculated in pseudodiatomic model assuming fixed values 

length and HNH angle, R = 1.987 a , 9 = 101.14°. 

1 o 


of NH bond 
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In the quasi-dlatomic model used here, the 2 2 A" state is quasi-bound, 
supporting two levels, and the photodissociation cross section shows very nar- 
row resonance peaks centered at 1680 A and 1770 A with halfwidthtf of several 
cm -1 . However, additional MCSCF calculations show that the barrier to photo- 
dissociation disappears as the bond angle increases. In a one-dimensional 
model, a better estimate for the 2 2 A" cross section results from an 
artificial barrierless potential, assuming dissociation takes place along an 
effective reaction coordinate. Using such an artificial potential and the 
calculated transition moment function results in an estimated cross section 
for the 2 2 A" state with a maximum of ~ 1 x 10 ^ cm 2 at 1780 A and a 
halfwidth of 130 A. A more sophisticated treatment of photodissociation ions 
is required to describe the detailed structure of the cross section. 

The Importance of these photodissociation cross sections in modeling 
cometary atmospheres depends on both the magnitude of the cross sections and 
the transition wavelength at which they peak. In the wavelength region where 
NH 2 transitions have their maxima, the solar flux varies a great deal, 
increasing by four orders of magnitude as wavelength is increased from 1500 A 
to 3000 A. 13 While the 2 2 A f state has the smallest cross section of the 
transitions reported here, it takes place at the longest wavelength of the 
transitions studied, where the solar flux is higher. This makes its role more 
nearly comparable to that of the other states than would be concluded from the 
cross sections alone. 

The cross sections reported here are comparable in magnitude to processes 
that are at present considered in cometary models^ ^ and thus should reason- 
ably be included in future modeling studies. 

Calculations on C 2 H 

17 1 Q 

Based on the spectroscopy of the C 2 molecule, ' thermochemical data 10 for 
the dissociation energy of C 2 H, and previous calculations for C 2 H in the 
equilibrium geometry,^ we have drawn the correlation diagram shown in 
Figure 1. We are interested in studying dissociation of C 2 H to states of C 2 
up to the d 3 II states, which is the upper state for Swan bands. We are 

O 

considering only doublet states of C 2 H, i.e., those that are dipole connected 
to the C 2 H X 2 Z + ground state. Therefore, thi6 work is concentrated on the 
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first two C 2 H states of ^£ + symmetry, the first state of symmetry, and the 

o 

’first three states of It symmetry. All of these states are valence in 

2 

character except the 3 It state, which appears to be a Rydberg state. 

Extensive MCSCF calculations have been performed for the five valence 

states of interest. A basis set appropriate for describing the Rydberg states 

as well as the valence states has been used from the outset. MCSCF potential 

curves for these five valence states in linear geometry are shown in 

Figure 2. For geometries near the ground state equilibrium geometry, the 

dominant configurations for the 2 £ , 2 H, and 1 A £ states all have the ir 

orbital occupied. For geometries where C 2 is separated from the H atom, the 
4c 

it orbital is not occupied. Thus, the potential curves in Figure 2 for these 
states Involve curve crossings that result in the barrier to linear dissocia- 
tion shown here. Therefore, we have undertaken calculations for bent geome- 
tries to determine the lowest energy path to dissociation. Photodissociation 
would then be treated using an effective reaction coordinate. 


JA-330583-275 


FIGURE 2 


POTENTIAL ENERGY IN ATOMIC UNITS FOR LINEAR 
C 2 H FOR FIXED C-C DISTANCE - 2.304 a 0 AS A 
FUNCTION OF C-H BOND LENGTH. 
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